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ABSTRACT 

Estimations of black hole spin in the three Galactic microquasars GRS 1915+105, GRO J1655-40, and XTE J1550-564 have been 
carried out based on spectral and timing X-ray measurements and various theoretical concepts. Among others, a non-linear resonance 
between axisymmetric epicyclic oscillation modes of an accretion disc around a Kerr black hole has been considered as a model for 
the observed high-frequency quasi-periodic oscillations (HF QPOs). Estimates of spin predicted by this model have been derived 
based on the geodesic approximation of the accreted fluid motion. Here we assume accretion flow described by the model of a 
pressure-supported torus and carry out related corrections to the mass-spin estimates. We And that for dimensionless black hole 
spin a = cJ/GM 2 < 0.9, the resonant eigenfrequencies are very close to those calculated for the geodesic motion. Their values 
slightly grow with increasing torus thickness. These findings agree well with results of a previous study carried out in the pseudo- 
Newtonian approximation. The situation becomes different for a > 0.9, in which case the resonant eigenfrequencies rapidly decrease 
as the torus thickness increases. We conclude that the assumed non-geodesic effects shift the lower limit of the spin, implied for 
the three microquasars by the epicyclic model and independently measured masses, from a ~ 0.7 to a ~ 0.6. Their consideration 
furthermore confirms compatibility of the model with the rapid spin of GRS 1915+105 and provides highly testable predictions of 
the QPO frequencies. Individual sources with a moderate spin (a < 0.9) should exhibit a smaller spread of the measured 3:2 QPO 
frequencies than sources with a near-extreme spin (a ~ 1). This should be further examined using the large amount of high-resolution 
data expected to become available with the next generation of X-ray instruments, such as the proposed Large Observatory for X-ray 
Timing (LOFT). 
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1. Introduction 

Studying the X-ray spectra and variability provides a useful tool 
for probing strong-field gravity effects and places constraints on 
the properties of compact objects, such as the mass or spin of 
a black hole. One of the commonly accepted ways to measure 
the black hole spin is related to fitting the X-ray spectral contin¬ 
uum or the relativistic ally broadened Fe K alpha lines. Using the 
spectral fitting methods, estimations of black hole spin have been 
carried out by many authors in the past (see e.g., McClintock et 
al„ 2006,2008; Middleton et al„ 2006; Done et al„ 2007; Miller, 
2007; Shafee et al„ 2008; McClintock et al„ 2010, 2011, 2014). 

There is another approach to estimating parameters of ac¬ 
creting black holes (BH) and neutron stars (NS) that has been 
gaining popularity in the past decade. It is the determination of 
compact object properties through the theory of high-frequency 
quasi-periodic oscillations (HF QPOs). These peaked features 
have been observed in the X-ray power density spectra of the 
BH and NS low-mass X-ray binaries for several decades (see 
e.g., van der Klis, 2006; Belloni & Stella, 2014, for a review). In 
the BH systems, the HF QPOs appear at frequencies that often 
form rational ratios with a preferred ratio of 3:2 (Abramowicz & 
Kluzniak, 2001; McClintock & Remillard, 2003). 


1.1. Orbital models of HF QPOs 

A rich variety of models has been proposed to explain the 
BH HF QPOs. Many of them are associated with motion of 
matter accreting onto the central compact object (e.g., Kato 
& Fukue, 1980; Alpar & Shaham, 1985; Lamb et al., 1985; 
Okazaki, Kato & Fukue, 1987; Nowak & Wagoner, 1991; Perez 
et al., 1997; Kluzniak, 1998; Miller, Lamb & Psaltis , 1998; 
Psaltis et al., 1999; Stella & Vietri, 1998a,b; Wagoner, 1999; 
Wagoner et al., 2001; Silbergleit, Wagoner & Ortega-Rodrfguez, 
2001; Abramowicz & Kluzniak, 2001; Titarchuk & Wood, 2002; 
Rezzolla, Yoshida & Zanotti , 2003; Petri, 2005; Zhang, 2005; 
Sramkova et al., 2007; Kato, 2007; Cadez et al., 2008; Stuchlfk 
et al., 2008; Horak et al., 2009; Mukhopadhyay, 2009; Lai et al., 
2012; Ortega et al., 2014, and others). Despite the different in¬ 
volved physical mechanisms, several of the HF QPO models as¬ 
sign the observed frequencies to orbital frequencies that can un¬ 
der some approximations be expressed in terms of Keplerian and 
epicyclic frequencies of perturbed circular geodesic motion. 

One of the early models, for instance, the so-called 
relativistic-precession model (Stella & Vietri, 1998a,b, 1999), 
identifies the kHz QPO frequencies v v and v, with the Keplerian 
orbital frequency and the periastron-precession frequency, v v = 
v K , v L = v K - v r (where v r denotes the radial epicyclic frequency). 
Within the model framework, it is usually assumed that the ob¬ 
served kHz variability originates in a bright localized spot or 
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Table 1 . Properties of the three microquasars GRO 1655-40, GRS 1915+105, and XTE 1550-564. The individual columns display the frequencies 
of the lower and upper QPO peaks (Strohmayer, 2001a,b; Remillard et al., 2002, 2003), the mass estimates (Greene et al., 2001; Greiner et al., 
2001; Orosz et al., 2002; McClintock & Remillard, 2003, 2006), and the spin values inferred from the epicyclic resonance model assuming a 
geodesic accretion flow. 


Source 

vl [Hz] 

vu [Hz] 

Mass [Me] 

a 

GRO 1655-40 

300 

450 

6.0-6.6“ 

0.96-0.99 

GRS 1915+105 

113 

168 

10.0-18.0* 

0.68-0.99 

XTE 1550-564 

184 

276 

8.4-10.8 

0.89-0.99 


Notes. (a) Alternative predictions of the black hole mass in GRO 1655-40 are 5.1 - 5.7 M 0 (Beer & Podsiadlowski, 2002) and 5.8 - 6.8 M 0 (Greene 
et al., 2001). {b) For GRS 1915+105, Reid et al. (2014) have recently obtained a tighter mass estimate, 10.6 - 14.4 M G , based on a trigonometric 
parallax measurement of the distance (see also Steeghs et al., 2013). 


blob orbiting the black hole at a slightly eccentric orbit. The pho¬ 
tons emitted by the spot, which have a high orbital velocity of the 
order of percents of the speed of light, are then propagated from 
the strong gravity region towards the observer. Consequently, the 
detected X-ray flux is periodically modulated by the relativistic 
effects (see e.g., Schnittman & Bertschinger, 2004; Bakala et 
al., 2014, and references therein). 

A different class of models is based on the assumption that 
the light curve reflects variations in the emission of an accretion 
disk. A theory of normal modes of accretion disk oscillation, 
referred to as discoseismology, has been developed in the pa¬ 
pers coauthored by Kato and by Wagoner. For instance, Kato & 
Fukue (1980) consider the g-mode of disk oscillation, whose fre¬ 
quency is close to the highest radial epicyclic frequency, while 
Silbergleit, Wagoner & Ortega-Rodrfguez (2001) developed a 
theory of c-mode (disc corrugation) whose frequency is related 
to the difference between the vertical epicyclic frequency and the 
orbital frequency. 

1.1.1. Resonance models 

Kluzniak & Abramowicz (2001) have introduced another con¬ 
cept of orbital models. In these models, the relativistic effects 
acting on the emitted photons are also essential, but the HF 
QPOs are attributed to a non-linear resonance between two 
modes of accretion disc oscillations. This well explains the ob¬ 
served preferred frequency ratio, Vu/v L = 3/2. Since the concept 
deals with a collective motion of accreted matter, it has a high 
potential of the observed flux modulation (Bursa et al., 2004; 
Mazur et al., 2013; Vincent et al., 2014; Bakala et al., 2014). 

The framework of resonance models has been extensively 
developed (Abramowicz & Kluzniak, 2001; Abramowicz et al., 
2003; Kluzniak et al., 2004; Horak & Karas, 2006; Horak et al., 
2009) and miscellaneous combinations of disc oscillation modes 
have been considered. In several particular cases, the observed 
frequencies can be expressed in terms of Keplerian and epicyclic 
frequencies (see Torok et al., 2005; Horak, 2008; Kluzniak, 
2008, for details and references). 

1.2. Spin estimation and epicyclic resonance model 

In Kerr geometry, which describes the spacetime around a rotat¬ 
ing black hole, the frequencies of the radial and vertical epicyclic 
oscillations depend on the mass and spin of the black hole. It is 
therefore possible to determine the BH mass or spin from the 
epicyclic resonance model and the observed 3:2 QPO frequen¬ 
cies. For three Galactic microquasars with the 3:2 HF QPOs 
(GRS 1915+105, GRO J1655-40, and XTE J1550-564), such 


BH spin estimations have been carried out in the past for the 
epicyclic resonance model and for several other HF QPO mod¬ 
els (Strohmayer, 2001a,b; Wagoner et al., 2001; Abramowicz 
& Kluzniak, 2001; Kato, 2004a,b; Torok et al., 2005, 2011; 
Motta et al., 2014a,b). Analogous spin estimations have been ob¬ 
tained for the three sources through other methods by different 
groups of authors using the X-ray continuum-fitting procedure 
(McClintock et al., 2006, 2008; Middleton et al., 2006; Miller et 
al., 2009), and the iron line method (Miller, 2007; Davis et al., 
2006; McClintock & Remillard, 2009; Blum et al., 2009; Steiner 
et al., 2011). The spin estimates obtained through the spectral 
and timing approaches are as yet somewhat inconsistent, how¬ 
ever (see a detailed discussion in Torok et al., 2011). 

In this work, we consider one of the popular versions of res¬ 
onance models represented by the epicyclic parametric (or inter¬ 
nal) resonance model that assumes a 3:2 non-linear resonance 
between the axisymmetric radial and vertical epicyclic modes of 
accretion disc oscillations. All the estimates derived so far that 
were related to the epicyclic resonance model assume a geodesic 
approximation of the accreted fluid motion. Within this approxi¬ 
mation, the two observable resonant frequencies are represented 
by the radial and vertical epicyclic frequencies, v r and vg, of test 
particle motion. In more general accretion flows, non-geodesic 
effects connected for instance to pressure gradients, magnetic 
fields, or other forces may affect the properties of the consid¬ 
ered oscillation modes and consequently the inferred spin pre¬ 
dictions. 

In the following, we explore black hole spin predictions 
implied by the epicyclic resonance QPO model consider¬ 
ing non-geodesic influence that originates in pressure forces 
present in accretion flow, which is modelled by a small 
equilibrium-pressure-supported, perfect-fluid torus. Properties 
of the epicyclic modes of torus oscillations, such as modifica¬ 
tions to their frequencies due to pressure gradients present in 
the torus, were calculated by Blaes et al. (2007) in the pseudo- 
Newtonian approximation and were later generalised by Straub 
& Sramkova (2009) for Kerr geometry. Using the results of 
Straub & Sramkova (2009), we apply the modified epicyclic fre¬ 
quencies to calculate the corresponding estimates of black hole 
spin based on the epicyclic resonance QPO model. 

2. Resonant eigenfrequencies and black hole 
parameters implied by the epicyclic model 

The explicit formulae of the epicyclic frequencies in Kerr ge¬ 
ometry were first derived by Aliev & Galtsov (1981) and may 
be written in the Boyer-Lindquist coordinates f, r. 0, <p as (e.g., 
Silbergleit, Wagoner & Ortega-Rodrfguez, 2001; Torok and 


2 








E. Sramkova et al.: Black hole spin inferred from 3:2 resonance model 


Stuchlfk, 2005) 


2 2 
V r = arVjr, 

Vg = a e v 
where 


v K = 


GM 


G 2 


1/2 


(x 3/2 + fl) 


-1 


a r (x ,a) 5 1 — 6 x~' + 8 a x” 3 ^ 2 -3 a 2 x~ 2 , 
ag (x, a) = 1 - 4 a x -3 / 2 + 3fl 2 x~ 2 , 
x = r/rc, rG = GM/c 2 . 


(1) 

( 2 ) 


(3) 

(4) 

(5) 

(6) 


Here M and a denote the mass and the dimensionless angular 
momentum (spin) of the black hole. 

Within the epicyclic resonance (ER) model, it is assumed 
that the resonant eigenfrequencies are equal to epicyclic frequen¬ 
cies defined at the xyi orbit where we have 


= 3/2v r . (7) 

Furthermore, it is assumed that the eigenfrequencies are equal 
to the observed upper and lower 3:2 QPO frequencies, v„ and 
Vi (e.g., Torok et al., 2005). Using this set of formulae and the 
observed 3:2 HF QPO frequencies along with the independently 
estimated black hole masses, the intervals of spin predicted by 
the ER model were previously calculated for the three Galactic 
microquasars (GRS 1915+105, GRO J1655-40, and XTE J1550- 
564 - see Torok et al., 2005, 2011). The obtained results and rel¬ 
evant properties of the three sources are summarized in Table 1. 

2.1. Epicyclic frequencies in a slightly non-slender torus 

We assumed a non-geodesic flow that was modelled by an equi¬ 
librium, slightly non-slender, pressure-supported perfect-fluid 
torus. The torus was assumed to have a constant specific angular 
momentum distribution and orbit a rotating Kerr black hole. We 
assumed that the torus is formed by a perfect fluid that fills up 
closed equipotential surfaces. One of these surfaces has a cusp 
in the equatorial plane, which embodies features similar to those 
of the LI Lagrange point. The topology of the equipotential sur¬ 
faces and the equilibrium torus is illustrated in Fig. I. 1 . 

In the investigated non-geodesic flow, the radial and verti¬ 
cal epicyclic oscillations of the fluid are modified by pressure 
forces. Using a perturbation theory, within the accuracy of a 
second order in the thickness of the torus, Straub & Sramkova 
(2009) calculated explicit formulae for the pressure corrections 
to epicyclic frequencies in a slightly non-slender, constant spe¬ 
cific angular momentum torus orbiting a Kerr black hole. The 
calculated epicyclic frequencies can be written in the following 
form: 

v* = v r + p 2 C r (M, r c , a), C r (M , r c , a) < 0, (8) 

Vg = yg + p 2 Cg(M, r c , a), C g (M , r c , a) < 0, (9) 

where C r (M, r c , a) and Cg(M, r c , a) denote the negative pressure 
corrections evaluated at the centre of the torus, r = r c . Explicit 
forms of expressions for the epicyclic frequencies (8) and (9) are 
given in formulae (52) and (56) of Straub & Sramkova (2009). 

1 For a detailed description of the model of the torus and further ref¬ 

erences see Straub & Sramkova (2009). 



rsin(0) 


Fig. 1. Illustration of equipotential surfaces of the accretion torus. The 
dotted horizontal line denotes the equatorial plane. The yellow region 
illustrates tori of various thickness. The orange region, together with the 
yellow, illustrates a torus filling the equipotential surface with a cusp. 


The dimensionless parameter ft determines the torus thick¬ 
ness. It is defined as (Blaes, 1985) 




2 net 


r 2 2 jj2 ’ 
' c v K.c u 


( 10 ) 


where c s , v Kc . and U are the sound speed, Keplerian frequency, 
and /-component of the contravariant four-velocity of the perfect 
fluid, all defined at the centre of the torus, r = r c . Within the per¬ 
turbative approach applied by Straub & Sramkova (2009), it is 
expected that 0 < ft 1 «c 1. For the critical p values corresponding 
to tori with a cusp, the limit of/I —» 0 corresponds to an infinitely 
slender torus located at the marginally stable circular orbit. High 
values of p < 1 then correspond to very large tori whose inner 
edge approaches the location of the marginally bound circular 
orbit. The behaviour of epicyclic frequencies for tori of differ¬ 
ent thickness,/! € {0, 0.1, 0.2, 0.4), and a Kerr black hole with 
a e {0.0, 0.5, 0.9, 0.99) is illustrated in Fig. 2. 


2.2. Behaviour of commensurable frequencies for 
moderately and rapidly rotating Kerr black holes 

Using formulae (8) and (9), we investigated the behaviour of 
epicyclic frequencies in fluid tori when they are in the 3:2 ratio. 
For Schwarzschild and moderately rotating Kerr black holes, the 
resonant frequencies behave qualitatively in the same fashion. 
The resonant radius in this case decreases with growing torus 
thickness, while the individual epicyclic frequencies defined at 
this radius increase with growing torus size. This is illustrated 
in Fig. 3, which, for different values of a and p , displays the 
ratio R of the vertical to radial epicyclic frequency, and the verti¬ 
cal epicyclic frequency calculated at the resonant radius. Panels 
a) and c) of Fig. 3 drawn for a — 0 and a = 0.5 show that 
from certain value of /?, there suddenly appear two different tori 
whose epicyclic frequencies display the 3:2 commensurability. 
This behaviour persists as we continue to increase the torus size 
until at some point the 3:2 resonant radius no longer exists. This 
phenomenon manifests itself in the form of a small loop in the 
profile of the resonant vertical epicyclic frequency plotted over 
P (panels b) and d) of Fig. 3). The same behaviour was found 
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a) b) 




Fig. 2. Epicyclic frequencies calculated at the centre of the torus, r = r c , plotted for various torus thicknesses (J3) and black hole spins a. The curves 
are drawn solely within the region where the considered value of /? does not exceed the critical value corresponding to a torus with a cusp. 


by Blaes et al. (2007), who carried out a similar analysis in the 
pseudo-Newtonian approximation. We find that for Kerr black 
holes, two different values of resonant frequency are allowed for 
the same value of /3 when a < 0.86. 

For higher values of a, the behavioural nature of the resonant 
frequencies begins to change significantly. The resonant radius 
retains its decreasing pattern with growing torus size, but two 
different tori with the 3:2 commensurability are no longer pos¬ 
sible. Rather distinct modifications also appear in the character 
of the resonant frequency, which at first increases with rising /j, 
but at some point starts to decrease rapidly. For a black hole spin 
approaching a = 0.89, the resonant frequencies implied by the 
highest allowed values of yS become lower than those calculated 
for the geodesic motion. This is illustrated in panels e) and f) of 
Fig. 3. For a still higher black hole spin (a > 0.99), characteris¬ 
tic features of the resonant frequencies continue to change even 
further. While the resonant radius continues to drop monotoni- 
cally with growing torus thickness, resonant frequencies are no 
longer higher than those of free test particles. They now decrease 
for tori of all sizes (J3 > 0), see panels g) and h) of Fig. 3. This is 
because the nature of epicyclic frequencies for rapidly rotating 
black holes is such that the decrease of frequencies due to higher 
/? overrides the increase implied by lower resonant radius. 

2.3. Comparison of predicted 3:2 epicyclic frequencies with 
observation 

Figure 4 uses the relativistic 1 /M scaling of the orbital frequen¬ 
cies. It displays the Mxv v (a ) relation expected from the geodesic 
epicyclic resonance model marked by the thick black line. The 
shaded region describes the mass-spin interval predicted by the 
epicyclic resonance model for a slightly non-slender torus of 
various thickness. In the figure, we also present the conserva¬ 
tive observational bounds on M x v v determined for the three 
microquasars (see Table 1). These bounds are indicated by the 
colour-coded areas. Apparently, while the geodesic line enters 
the observationally determined vertical range of microquasars 
when a > 0.7, the non-geodesic line enters this range already for 
a > 0.6. 2 


2 According to individual studies based on spectral and timing meth¬ 
ods, it would also be possible to include several (vertical) bounds on the 
spin of microquasars. Within the figure, however, we focus primarily on 
the epicyclic resonance model of HF QPOs and its predictions. The dis- 


Figure 4 clearly illustrates that for a < 0.9, the highest pos¬ 
sible increase of the predicted resonant frequency due to pres¬ 
sure effects is rather small, reaching about 10%. However, for 
a > 0.9, the shaded region in the figure shifts from the area 
above the /3 - 0 curve to the area located mostly below this 
curve. The strongest possible shift of the resonant frequency is 
now much greater than it was for a < 0.9. The coloured curves in 
Fig. 4 illustrate the mass-spin values predicted by the model for 
tori of p — 0.3, /? = 0.4, and /3 = 0.5. In this case, predictions of 
the model with non-geodesic flow may significantly differ from 
those calculated for the geodesic flow. This has a direct link to 
observations, since the calculated shaded region overlaps a large 
part of the observationally determined yellow area. 


3. Discussion and conclusions 

Following the previous studies of Blaes et al. (2007), Straub & 
Sramkova (2009), and Tdrok et al. (2011), we explored here the 
implications of consideration of pressure forces that originate in 
a pressure-supported fluid torus for the predictions of black hole 
spin inferred from the epicyclic resonance QPO model. Our find¬ 
ings for moderately rotating Kerr black holes with a < 0.9 agree 
well with those carried out in the pseudo-Newtonian study of 
Blaes et al. (2007). The epicyclic resonance model with a non¬ 
geodesic flow predicts the black hole spin to be slightly lower 
than the spin estimated previously by the same model with a 
geodesic flow. The lower limit on the spin of Galactic micro¬ 
quasars inferred by the model is then decreased to a ~ 0.6. These 
findings indicate that for up to a ~ 0.9, the model has a high pre¬ 
dictive power (falsifiability). 

Different but very interesting results appear for rapidly ro¬ 
tating black holes with a > 0.9. In this case, the behaviour of 
the resonant frequencies begins to markedly change in compar¬ 
ison to the pseudo-Newtonian case. This occurs in such a way 
that the predicted frequencies can be both higher and lower than 
those predicted by the geodesic model, with the strongest pos¬ 
sible shift being much more significant towards the lower val¬ 
ues. Therefore, when the source mass is determined from the 
3:2 QPO frequencies, the epicyclic resonance model has a lower 


cussion of various spectral and timing estimates of a can be found for 
instance in Miller (2007), McClintock et al. (2011, 2014), Steiner et al. 
(2013), Torok et al. (2011), and Kotrlova (2014). 
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a) 
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c) 


d) 



e) 
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Fig. 3. Ratio R of the vertical to radial epicyclic frequency and the resonant vertical epicyclic frequency Vg calculated at the radius X 31 . The arrows 
in panels a)-d) indicate the behaviour corresponding to the two resonant radii jc 3:2 . The inner (red) and outer (blue) resonant radius is denoted. 


predictive power for rapidly rotating black holes than for moder¬ 
ately rotating black holes. On the other hand, our findings con¬ 
firm the compatibility of the model with the rapid spin of the 
microquasar GRS 1915+105 (McClintock et al., 2006; Torok et 
al., 2011; Fragos & McClintock, 2014), since the expected mass 
interval almost overlaps with the interval constrained observa- 
tionally. Furthermore, we provided highly testable predictions of 
the QPO frequencies. When the torus size (J3 ) is not determined 
by some additional physical requirements, a wide range of com¬ 
mensurable frequencies is predicted for near-extreme rotating 
black holes. 3 Individual sources with a moderate spin (a < 0.9) 
should therefore exhibit a smaller spread of the measured 3:2 
QPO frequencies than sources with a near-extreme spin (a ~ 1). 
At present, there is a rather small amount of data needed for 
further analysis of this issue. Our predictions, however, can be 


3 Values of /3 as high as /? = 0.5 considered in panel h) of Fig. 3 may 
require additional treatment, since these values exceed the limits of full 
applicability of the adopted approximation. Nevertheless, the illustrated 
phenomenon is very significant and well apparent already for ~ 0.3. 


explored using the large amount of high time-resolution data ex¬ 
pected to become available with the next generation of X-ray 
instruments, such as the proposed Large Observatory for X-ray 
Timing (LOFT; Feroci et al., 2012). 
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